Theory of tunnelling into and from cuprates 



A.S. Alexandrov 



Department of Physics, Loughborough University, Loughborough LE11 3TU, U.K. 

Abstract 

A single-particle spectral density is proposed for cuprates taking into ac- 
count the bipolaron formation, realistic band structure, thermal fluctuations 
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and disorder. Tunnelling and photoemission (PES) spectra are described, in 
eluding the temperature independent gap observed both in the superconduct- 
ing and normal states, the emission/injection asymmetry, the finite zero-bias 
conductance, the spectral shape in the gap region and its temperature and 

O . 

doping dependence, dip-hump incoherent asymmetric features at high voltage 

(tunnelling) and large binding energy (PES). 
'■ PACS numbers:74.20.-z,74.65.+n,74.60.Mj 
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The strong-coupling extension of the BCS theory based on the 1/A multi-polaron per- 
turbation technique firmly predicts the transition to a charged Bose liquid in the crossover 
region of intermediate values of the BCS coupling constant A There is a fundamental 
difference between the bipolaron theory of high-T c cuprates || and other theories involving 
real-space pairs (bosons) tightly bound by a field of a pure electronic origin. As empha- 
sised by Emery et al [[| such 'electronic' theories are a priori implausible due to the strong 
short-range Coulomb repulsion between two carriers. The direct (density-density) repulsion 
is usually much stronger than any exchange interaction. On the other hand, the Frohlich 
electron-phonon interaction can provide mobile intersite bipolarons in the Cu02 plane con- 
densing at high T c Q. The (bi)polaronic nature of carriers in cuprates explains a very 
small coherence volume in the superconducting state, the mid-infrared conductivity 0, the 
isotope effect on the carrier mass |J. 

Although the charged Bose liquid of bipolarons describes anomalous thermodynamics 
and kinetics of cuprates [J/J, finite frequency /momentum response functions of bipolaronic 
superconductors remain to be established. In this letter we derive a single-particle spectral 
function of strongly coupled carriers in a random potential which provides a quantitative 
description of recent tunnelling spectra and explains some photoemission features (see 
[T0|1 and references in |^-|10|). 

In the framework of our theory the ground state of cuprates is a charged Bose-liquid of 
intersite bipolarons with single polarons existing only as excitations with the energy A/2 
or larger. Different from the BCS description the pair binding energy A is temperature 
independent. Hence, there is no other phase transition except a superfluid one at T = 
T c . The characteristic temperature T* of the normal phase is a crossover temperature of 
the order of A/2 where the population of the upper polaronic band becomes comparable 
with the bipolaron density. Along this line the theory of tunnelling in the bipolaronic 
superconductors was developed both for two-particle ]TT[ and single- particle |L2] transitions 



through a dielectric contact. It allowed us to understand the temperature independent gap 
and the asymmetry of the current-voltage characteristics, observed already in the earlier 



tunnelling experiments However, an attempt to fit the conductance structure led us 

to a very narrow (bi)polaronic band with a bandwidth of the order or even less than T c 



TT| , [T2|| . Such a bandwidth is not compatible with the experimental estimate of the effective 
carrier mass, m* ~ 2 — 10m e (depending on doping) from the London penetration depth 0). 
It is also incompatible with the theoretical estimate [Q of the (bi)polaron bandwidth (c.a. 
lOOmeV or larger) based on the well-established value of the Frohlich interaction. Moreover, 
any description based on the Bloch representation is hardly justifiable for cuprates with the 
mean free path often comparable with the lattice constant. One has to consider a random 
potential and thermal fluctuations along with a strong pairing potential and band-structure 
effects. 

We apply a single-particle tunnelling Hamiltonian describing the injection of an electron 
into a single hole polaronic state P with the matrix element P^ v and into a paired hole 
(bipolaronic) state B with the matrix element B^^ (Fig.l), 

Htun = p ^ a ^Pv + E B^^akptbf, + H.c. (1) 

Here k, v and \x are the quantum numbers describing the electron in a tip (Fig.l) (annihila- 
tion operator at) a hole polaron (p„) and a hole bipolaron (b^) in the Cu02 plane in a random 
field, respectively. If the eigenstates \v) and \y) are known, the matrix elements P^ v and 
-Bk.z, lA t are derived by the use of the site representation and the canonical polaronic transfor- 
mation as discussed in detail in Ref. []IT , |12"| |. They are almost independent of k ~ kp, fi, v in 



a wide voltage and binding energy range, P^ v ~ P = constant, N 1 / 2 B^ V>11 ~ B — constant, 
N is the number of cells in the sample volume. In general, B and P are different |T2[ because 



the second hole in a small coherence volume changes the potential barrier of the contact for 
the tunnelling B compared with P. 

The injection rate is given by the Fermi Golden Rule as 

W in = 2nP 2 J2 + eV + + A/2) 

k,K 

+ ^ E /(&K[1 " f(^)Mu + eV + tr-&- A/2), (2) 



and the emission rate is 



W em = 2vrP 2 ]T[1 ~ /(&)][! - /(&)]<K£k + eV + &, + A/2) 

+ ^ E [1 - /(&)] K + l]/(60*(& + eV + 6. - fc, - A/2). (3) 

Here /(£) = [exp(£/T) + is the Fermi-Dirac distribution function for electrons in the 
tip and hole polarons in the sample, n M is the bipolaron distribution jTJJ , ^k,^ is the energy- 
spectrum of electrons, polarons and bipolarons, respectively, V is the sample voltage, e is 
the elementary (positive) charge, and h — ks — 1- Taking the difference of the injection 
and emission rates multiplied by e and differentiating it with respect to V we obtain the 
conductance a = dl/dV as 

7iN iin P 2 r°° o If + eV + A/2' 



a(V,T) = ^f- T d^p(Osech 2 
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+ J d£p(0 2>„ + /(Ol^ec/i 2 [e - e M - + A/2 o^er2T] , (4) 



with N tip and p(f ) the density of states (DOS) per spin in the tip near the Fermi level and the 
polaronic DOS in the sample, respectively. To arrive at the analytical result we assume that 
the Coulomb bipolaron-bipolaron repulsion [|7| is relatively weak, and temperature is small 
compared with A and the characteristic energy scale €q of the polaronic DOS (see below). 
Then the bipolaron distribution n M is narrow [1TJJ] and the population of the polaronic states 
is low. That allows us to integrate out the bipolaronic states taking £ M = /(£) = in the 
second term of Eq.(4). The derivative of the Fermi distribution, sech 2 (x) cuts the integrals 
at \x\ — 1. As a result we obtain 

a(V, T) = ^ p [N(-eV - A/2 + 2T) - N(-eV - A/2 - 2T)\ 
ite 2 xN+- B 2 

+ Up [N{eV - A/2 + 2T) - N{eV - A/2 - 2T)} , (5) 

with N(E) = Jl^ dE'p(E') the cumulative polaronic DOS, and x the bipolaron density per 
cell proportional to the doping. The first term in Eq.(5) proportional to P 2 contributes 
to the emission, while the second one describes the injection. They have different values 



because, in general, we have B > P and x < 1, so that A = P 2 /B 2 x differs from unity [12 . 



We notice that the spectral shape of the emission and the injection might be also different 
if the bipolaron distribution is sufficiently wide. There is an additional integration in the 
second term in Eq.(4), describing the injection compared with the emission. Hence, a fine 
spectral structure like the dip found in the emission ]T0|,P[ (and references therein) can be 
nearly washed out from the injection as observed 

The tunnelling (or PES) spectrum, Eq.(5) is defined if the polaronic cumulative DOS is 
known. It depends on the band structure, dressing and scattering. While high frequency 
phonons and magnetic fluctuations are responsible for the high-energy spectral features in 
the region of the order of the Franck-Condon (polaronic) level shift (> 100 meV) J7j, the low 
energy spectral function is determined by the band structure and by the thermal lattice, spin 
and random fluctuations. The p-hole polaron in cuprates is almost one-dimensional due to 
the large difference in the ppa and ppn hopping integrals and the effective 'one-dimensional' 
localisation by the random potential as described in Ref. 0]. This is confirmed by the 
angle-resolved PES [17 with no dispersion along certain directions of the two-dimensional 



Brillouin zone. Because the amount of disorder is high and the screening radius is about 
the lattice constant, we can describe the effect of disorder and of the thermal fluctuations as 
'white Gaussian noise'. The relevant spectral density A(k, E) for a one- dimensional particle 
in a random Gaussian potential was derived by Halperin [|15|] and the density of states, 
p{E) by Frish and Lloyd [[Uj. Halperin found the spectral function both for a 'Schrodinger' 
particle (i.e. in the effective mass approximation) and for a 'discrete' particle (tight-binding 
approximation). The estimated polaronic bandwidth is about 100 meV or larger, so we 
consider the effective mass (m*) approximation for the spectral density, given by |T5[ 
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with Po,i(z) obeying two differential equations 



Po = 0, (7) 
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and 

d 2 d 



. „ +-f (z 2 + 2E) -z-ik 
dz z dz 



Pi+p = 0. (8) 



The boundary conditions are to be found in Ref . |H| . Here the energy E and the momentum 
k are measured in the units eo = (D 2 m*Y^ 3 and ko = (.D 1 / 2 ??!*) 2 / 3 , respectively The constant 
D describes the second moment of the Gaussian potential comprising thermal and random 
fluctuations as D = 2(V 2 T/M + 

n im v o), where Vq is the amplitude of the deformation 
potential, M is the elastic modulus, rii m is the impurity density, and vq is the coefficient of 
the 5- function impurity potential. Then the cumulative DOS 

/E roo 
dE' dkA(k,E') (9) 
-oo J — oo 

is expressed analytically in terms of the tabulated Airy functions Ai(x) and Bi(x) as 

N(E) = tt~ 2 [Ai 2 (-2E) + Bi 2 (-2E)] ^ . (10) 

Substituting Eq (10) into Eq.(5) we obtain the expression for the tunnelling and PES spectra 
in the voltage (energy) region, where the high-frequency phonon shake-off is forbidden by 
the energy conservation. In particular, we find for T = 

with 

4 Ai(~2x)Ai'(-2x) + Bi(-2x)Bi'(-2x) 
Py x > ~ ^2 x [Ai 2 {-2x) + Bi 2 (-2x)] 2 * ( ' 

and the constant do = 2iik eQ 1 e 2 xNu p B 2 . 

We compare the conductance, Eq.(ll) with the scanning tunnelling microscope (STM) 
|| and point-contact tunnelling (PCT) 0] measurements in an overdoped and optimally 
doped Bi 2 Sr 2 CaCu20s+5 (Bi-2212) in Fig. 2 and Fig. 3, respectively. The bipolaron theory 
describes quantitatively the spectra in the gap region, including the zero-bias conductance 
at T = 0, the asymmetry, and the decreasing background at higher voltages which are in- 
consistent with the classic BCS theory, no matter s or d-wave. The zero-bias conductance 



at T = is explained by the presence of the impurity tails of the polaronic DOS inside the 
gap, while the decreasing background, proportional to | | x / 2 for \V\ >> A/2 is explained 
by the one-dimensional band dispersion of polarons. The peak amplitudes and the zero-bias 
conductance are determined by the ratio of the bipolaron binding energy A to the charac- 
teristic scattering rate 6q. The position of two peaks and their amplitudes relative to the 
zero-bias conductance allow us to determine the relevant parameters A, eo and the asym- 
metry A with an error bar less than five percent. The doping dependence of A agrees with 
that found from the uniform magnetic susceptibility and explained by us as resulting 
from the screening of the Frohlich interaction by free carriers. 

An essential feature of the bipolaron theory is the temperature independent gap with the 
ratio A/T c , which might be quite different from the BCS one, 2Abcs/T c — 3 — 5. We obtain 
a very large ratio, A/T c ~ 9 for both overdoped, Fig.2 and optimally doped, Fig. 3, samples. 
Such a large ratio is difficult to understand in the framework of the BCS theory including 
its canonical strong-coupling extension. Renner et al || emphasised that the evolution 
of the STM spectra with temperature was very different from the classic BCS behaviour, 
both s and o?-wave. The 'superconducting' gap was found temperature independent evolving 
into the 'normal' gap above T c . In our theory this is one and the same gap, which is the 
bound energy of real-space bipolarons. It does not disappear at any temperature, neither 
at T c nor at T*. The theoretical evolution of the spectrum with temperature described in 
Eq.(5), Fig. 3 (inset), reproduces well the experimental features. In particular, the zero-bias 
conductance increases with temperature and there is no sign that the gap closes at a given 
temperature, in agreement with the experiment ||. The latter observation rules out any 
role of superconducting phase or spin fluctuations in the normal gap. We notice that the 
theoretical peaks shift to higher energies above T c , Fig. 3 (inset), as observed (see Fig.2 in 
Ref. 0). 

There is some characteristic voltage (binding energy) V c , (Fig.2 and Fig.3), above which 
the experimental STM and PCT conductance deviate from the theoretical one. We believe 
that a hump observed above V c is due to a polaronic cloud, as discussed in Ref. @]. The 
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high-frequency phonons and magnetic fluctuations contribute to the excess spectral weight 
with a maximum around twice the Franck-Condon shift. The dip of the spectral weight 
at V c is explained by the electron-collective excitation coupling as suggested by Shen and 
Schrieffer JTD] for PES and discussed by DeWilde et al || for STM and PCT. We have shown 
earlier that a similar dip structure appears in the electronic DOS as a result of the electron- 



Einstein phonon interaction [19||. With the established DOS, Eq.(12) we can quantify this 
feature. If we determine V c ~ 70meV in the overdoped sample, Fig. 2, and V c ~ 80meV in 
the optimally doped, then the relevant polaron kinetic energy, E c = eV c — A/2, appears to be 
about 39meV and 42meV, respectively. Hence, the emission of any dispersionless phonon of 
that energy can produce a dip in the polaronic DOS as described in Ref. [pl|. However, E c 



is also close to the estimated energy of the magnetic fluctuations thought to be responsible 
for the spin-flip neutron scattering at q — (jr, it) ( 4lmeV peak observed in Y BCI2CU3O7 
p0|). These magnetic fluctuations, should they be found in Bi2212, might contribute to the 
dip as well. The polaron-bipolaron inelastic collisions also contribute to the dip in a manner 



similar to the polaron-phonon interaction [TSj because V c is close to A. The polaronic DOS 



is finite in the middle of the gap, as discussed above. Hence, the polaron with the kinetic 
energy E c can break a bipolaron into two single polarons with the energy about A/2 each. 

The present theory of tunnelling and PES can be generalised to describe SIS junctions, c— 
axis current at high voltage and the angle- resolved photoemission (ARPES). In particular, 
SIS current- volt age characteristics are obtained by a convolution of the polaronic DOS with 
itself. As a result we get two peaks around the temperature independent \eV\ = A, and 
the symmetric shape of the SIS conductance, as observed ||. We expect an S'-shape / — 
V for the c — axis current both in the superconducting and normal state as a result of 
the bipolaron-breaking at high voltage. The S shape c-axis I — V characteristics were 
measured in Bi2212 and interpreted as a 'subgap resistance of intrinsic Josephson junctions' 
between superconducting Cu02 planes [^TJ. We believe, that while the notion of 'subgap' 
is perfectly correct, one cannot refer the observed nonlinearity to the Josephson junctions. 
The nonlinearity is still found at temperatures as high as 140K (Fig. 4 of Ref. |2l|) where 
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the C11O2 planes are definitely in the normal state, as expected in the framework of the 
bipolaron theory. Therefore, there is no doubt that the true normal state c — axis resistivity 
and the upper critical field H C 2 were measured in Ref. |22j by applying an external magnetic 
field. ARPES can be described with the spectral function A(k,E), determined numerically 
from Eqs (7,8). While such a feature of ARPES as the normal state gap is understood within 
the present analysis, the k dispersion will be presented elsewhere. 

In summary, we propose a single-particle spectral function for cuprates, which describes 
the spectral features observed in tunnelling and photoemission, in particular the temperature 
independent gap and the anomalous gap/T c ratio, injection/emission asymmetry both in the 
magnitude and in the shape, zero-bias conductance at zero temperature, and the spectral 
shape inside and outside the gap region. The dip- hump features and the temperature/doping 
dependence of the tunnelling conductance and PES are discussed as well. 

The author highly appreciates enlightening discussions with G.S. Boebinger, A.M. Camp- 
bell, R.A. Doyle, O. Fisher, W.Y. Liang, Ch. Renner, S.G. Rubin, J.R. Schrieffer, Z.-X. 
Shen, M. Springford, and V.N. Zavaritsky. We are grateful to M.A. Alexandrov for his 
assistance in computer calculations. We thank E.N. Sladkovskaia for her careful reading of 
the manuscript. 
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Figure Captures 

Fig.l SIN tunnelling into a single polaron state P and into a bipolaron state B. 

Fig. 2 Theoretical tunnelling conductance (line) compared with the STM conductance 
(T = 4.2K) in overdoped Bi2212 (T c = 74.3K §). 

Fig. 3 Theoretical tunnelling conductance (line) compared with the PCT conductance 
(T = 4.2K) in optimally doped Bi22\2 (T c = 95K 0). Inset shows the temperature 
dependence of the conductance calculated with a constant eo in the gap region. 
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